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Abstract 

We calculate expressions for the quadrupole moments of nonstrange baryons in which the number 
of QCD color charges is Nc- Using only the assumption of single-photon exchange, we obtain 4 
relations among the 6 moments, and show how all of them may be obtained from Qa+p to 
0(1/ N^) corrections. We compare to the Nc = 3 case, and obtain relations between the neutron 
charge radius and quadrupole moments. We also discuss prospects for the measurement of these 
moments. 
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I. INTRODUCTION 



The most interesting observable properties of hadrons, such as masses, electromagnetic 
moments, and scattering amplitudes, fall squarely in the nonperturbative regime of QCD 
and thus resist first-principles analytic calculation. While lattice simulations for static ob- 
servables continue to improve numerically, they do not typically address questions regarding 
which features of QCD explain which various aspects of hadronic observables. Other tech- 
niques, such as the operator product expansion, relate perturbative quantities, hadronic 
observables, and nonperturbative matrix elements pertaining to vacuum structure. The 
1/Nc expansion []T|, with A^^^ the number of color charges, is alone among QCD-inspired 
techniques in providing a perturbative parameter at all energy scales. While the notion of 
varying the fixed value of A^c = 3 may seem peculiar on first glance, the l/N^ expansion 
is made phenomenologically useful by the observation that the structure of the physical 
universe would not be qualitatively very different were the precise value of A^c = 5, 11, or 

113 i. 

The 1/Nc expansion has proved quite useful in describing the properties of baryons, such 
as ground-state masses 0, magnetic moments [|, |^, and axial couplings as well as 

excited baryon masses [0) H) 0! axial couplings [|, |10|, |Tl|, and pion-production [1^ and 
electromagnetic decay and production ||13[ amplitudes. Similar studies have investigated 
baryons containing a heavy quark |]T^ 0, and the nature of the nucleon-nucleon interac- 
tion |T^, as well variety of other properties. 

The analysis proceeds by studying interaction Hamiltonian operators of known spin-flavor 
transformation properties acting upon the baryon states. Each such operator is characterized 
by a coefficient proportional to a particular power of l/N^. The latter may be obtained 
by counting the minimum number of gluon exchanges necessary to connect the requisite 
number of quark lines in order for the operator to have a nonvanishing expectation value for 
the baryon state, since as oc 1/A^c- Of course, the possibility that contributions from a given 
operator may add coherently over the quarks, giving an additional enhancement of 0{Nc) 
to the full matrix element, must be taken into account. To this scheme one must add the 
observation that not all apparently different spin-flavor operators are linearly independent 
on a given multiplet of baryons. Often there appear relations between seemingly distinct 
operators due to the structure of the spin-flavor group [SU(6) for three quark flavors] or 
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due to the fact that one considers operators acting upon a particular representation of 
that group (the ground-state 56 for Nc = 3, for example). Application of such operator 
reduction rules 0, are then necessary to remove redundant degrees of freedom; see 
Sec. |ITT| for explicit examples. Together, these effects provide a power-counting scheme for 
deciding which Hamiltonian operators are most significant in determining contributions to 
the physical properties listed above. 

While we have been discussing the baryon as an A'^c-quark state, of course the physical 
baryon is a much more complicated object, containing a complex soup of gluons and sea 
quark-antiquark pairs. The simple-minded picture may be justified by noting that physical 
baryons for A^^ = 3 have precisely the same quantum numbers as those predicted by the 
simple quark model. Therefore, one can invent interpolating fields whose effect is to subsume 



the entire complicated substructure into Nc precisely defined constituent quark fields [18 



This identification allows the old constituent quark model to be placed on a rigorous footing 
and elucidates the nature of the physical picture used in the studies noted above. 

Once the relevant operators are ordered in a 1/Ac expansion and any operator coefficients 
due to known suppressions are taken into account [such as SU(3) flavor or isospin breaking], 
the matrix elements of the spin-flavor structures are computed by standard group-theoretical 
methods. Then, only an unknown dimensionless coefficient determined by strong interac- 
tion effects remains uncomputed . This approach is none other than the application of the 
Wigner-Eckart theorem in a novel context: The matrix elements of the spin-flavor structures 
constitute Clebsch-Gordan coefficients, and the leftover coefficients are reduced matrix ele- 
ments. Inasmuch as the expansion is natural, i.e., takes into account all known suppressions 
or enhancements, these reduced matrix elements should be of order unity. 

In the most general application of the method, one first generates a complete set of 
linearly independent operators, which is guaranteed to exist since the baryons belong to a 
finite representation for any particular Nc. Since ultimately one sets Nc 3, it is sufficient to 
consider only up to 3-body operators, i.e., irreducible operators involving 3 distinct quark 
lines. While this would seem to offer no predictivity — after all, then one would have an 
equal number of observables and operator coefficients — the hierarchy specified by the 1/Nc 
expansion determines which operators are more or less phenomenologically significant. 

In fact, the method can easily be adapted to various model assumptions. Suppose, 
for example, that one works in a model in which isospin violation arises from only one 
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source, say, quark charges. The most general expansion allows for every possible source of 
isospin violation; therefore, a number of possible operators can be neglected in this model. 
Since each independent operator is sensitive to a particular pattern of spin-flavor symmetry 
breaking that contributes only to certain combinations of baryon observables, the absence 
of each such operator in the model indicates a relation obeyed by the baryon states. This 
is precisely the same as the approach that was adopted for studying baryon charge radii in 
Ref. [|18], and for the Nc = 3 case for quadrupole moments in Ref. [1T9| . 

In this work we consider the quadrupole moments of the non-strange ground-state 
baryons. The quadrupole moments hold interest in providing vital information on the charge 
distributions of the baryons, which are inherently non-perturbative quantities. Only the non- 
strange baryons are considered in this work, as explained below, because the tricks of group 
theory for the 2-flavor case are different from those necessary to study the 3-flavor case. The 
latter case for quadrupole moments as well as charge radii will be considered in a subsequent 



But are baryon quadrupole moments measurable quantities? We address this central issue 



assumption is the one-photon exchange approximation, by which we mean that an arbitrary 
number of quarks can be involved in the quadrupole operator, but the photon itself couples 
to only one of them. The results of the analysis are presented in Sec. |V] and conclusions in 
Sec. |V[ 

II. PROSPECTS FOR MEASUREMENT 

The dearth of literature on baryon quadrupole moments compared to that on masses 
or magnetic moments points to the basic difficulty of measuring these quantities. The 
quadrupole tensor is rank-2, and therefore by the Wigner-Eckart theorem its expectation 
value, called the spectroscopic quadrupole moment Q, vanishes for the spin-1/2 nucleons. 
However, it should be noted [^] that Q, measured in space-fixed coordinates, is not identical 
to the second moment of the charge distribution. 



work |0|. 



in Sec. II. In Sec. ^TT| we generate the operator basis for a very general model: The only 




(2.1) 
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measured in body-fixed coordinates, with e the proton charge. Qo is called the intrinsic 
quadrupole moment, which would describe the shape of the charge distribution were it 
possible to take a snapshot of the particle at a fixed time. The most obvious manifestation 
of the intrinsic quadrupole moment of a spheroidal particle is whether it is prolate (football- 
shaped) or oblate (pancake-shaped). The relation between the two definitions for a particle 
of ground-state spin J is given by 



For J = or 1/2, the particle body-fixed axis points with equal probability in all directions 
by parity invariance, leading to the necessary vanishing of Q. However, Qq can be nonzero 
for such particles, and it has been shown in a variety of models 12^ that the proton is 



prolate and the is oblate. It has also been shown pSf that, given a sufficient number 
of measured electric quadrupole {E2) transition matrix elements, Qo can be extracted in a 
model-independent way. Nevertheless, the quadrupole moments discussed below are of the 
spectroscopic variety and are obtainable through quadrupole transitions. 

The spectroscopic quadrupole moments for the ground-state baryons consist of diagonal 
transitions for the spin-3/2 decuplet baryons A, S*, S*, and Q, and off-diagonal transitions 
between the spin-3/2 decuplet and spin-1/2 octet baryons, AA^, S*A, and Since 
all of the decuplet baryons except the fl decay strongly, couplings of the form decuplet- 
decuplet-7 can be measured only through virtual processes. Such experiments are difficult 



but not impossible; for example, Ref. p4| proposes the measurement of the A"*" magnetic mo- 
ment through the process 7p — > 'jir^p. However, the same authors note the near-impossibility 
of measuring the electric quadrupole moment through such experiments. The problem is 
that electric quadrupole {E2) operators are time-reversal odd. Thus, if the initial state 
and final state are identical (as for static diagonal quadrupole transitions), then this ma- 
trix element vanishes. There are only two ways around this restriction; the first is to use 
Coulomb photons rather than real radiation (the usefulness of Coulomb photons is limited 
by the short lifetime of the decuplet baryons), and the second is to extract the off-shell or 
recoil effects of the decuplet baryons, a proposition fraught with many difficulties. For our 
purposes, we suppose that the diagonal quadrupole transitions involving decuplet baryons 
will not be measured any time soon, and suggest that our calculation provides predictions 
of quantities that are very hard to measure. 
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The n~ is an exception to this constraint, since it has an appreciable hfetime. A number of 
clever experiments have been suggested to extract Qn, including the measurement of energy 
levels of an exotic atom formed by a heavy nucleus capturing an fi^ and the precession 
of spin as it traverses a crystal at small angles to the crystallographic plane PBI. In 
both cases, the purpose is to enhance the effect of the electric field gradient, to which the 
quadrupole moment couples. 

Since all the octet baryons (except the E") have appreciable lifetimes, the transition E2 
moments are all measurable in principle, inasmuch as one possesses a suitable source of the 
desired long-lived baryon. In the hyperon sector, the Primakoff reaction Y + Z Y* + Z , 
where Y{Y*) is the octet (decuplet) hyperon and Z is a heavy nucleus, is sensitive to both 



the Ml and E2 transition matrix elements; this process has been studied at SELEX ||2^ to 
obtain a bound on the radiative width of the S*. Experiments at Jefferson Lab that involve 
kaon photoproduction (7^ K^Y* K^Y'-f) ||28| can also provide useful information on 
these amplitudes. 

Lastly, the N —>■ A transition amplitudes have been studied in numerous experiments, 
for example see [2£, In particular, the transition quadrupole moment Qn^a may 

be extracted from the E2/M1 ratio (the relative quadrupole to dipole strength). A sam- 
ple recent measurement using pion photoproduction data is [—3.07 ± 0.26 (stat + syst) ± 
0.24(model)]% |2§, from which one obtains the value Qn->a = -0.108 ± 0.009 ± 0.034 
fm^. However, not all researchers agree on the size of uncertainties {e.g., those induced by 
a particular model used to describe the nonresonant background [^) when extracting the 
resonance amplitudes from the measurements. 



III. THE OPERATOR METHOD 



The operator method was described in the Introduction in general terms; here we lay 
out the specifics for this calculation. It is identical to the approach used for charge radii 
in |T^, and is the generalization (for nonstrange states) of the calculation performed in 
Ref. [l^. Any spin-flavor operator can be built from a basis using no more than three 1- 
body operators, i.e., spin-flavor operators acting upon one quark line. The most general such 
operators (restricted, of course, to those that transform as quadrupole tensors) are formed 
from sandwiching Pauli spin matrices and Gell-Mann flavor matrices (or isospin matrices in 
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the 2-flavor case) between quark creation and destruction operators. One then performs the 
operator reduction; for example, P = on the 2-flavor states. Were we to compute this 
model-independent expansion for the nonstrange baryon quadrupole moments, we would find 
precisely 6 operators, matching the number observables Qa++, Qa+, Qa", Qa-, Qa+p, and 
Qaou- Some of the operators in this list are suppressed by various powers of 1/Nc, leading 
to model-independent approximate relations between the quadrupole moments. Precisely 



this sort of expansion was carried out for nonstrange baryon charge radii in |118|. 

However, we adopt the mildly model-dependent but much more physical viewpoint that 
any quadrupole operator involves the quark charge only once. Operators containing, e.g., 
extra powers of the quark charge are suppressed by additional powers of e^/An = 1/137 
compared to the single-photon exchange approximation. This restricts the isospin structure 
of possible quadrupole operators considerably: Only one very particular combination of 
isosinglet and isovector coupling, given by the quark charge operator, then appears. 

The quadrupole spin coupling is a rank-2 tensor, and therefore clearly cannot be built 
with just one (rank-1) Pauli spin matrix. The diagonal rank-2 combination formed by two 
Pauli matrices acting on quarks labeled by i,j is just the familiar form: 

It is also possible to build a rank-2 spin tensor with 3 Pauli matrices, but such a form on any 
given quark pair can be divided into symmetric and antisymmetric Hermitian parts. The 
symmetric part turns out to be time- reversal odd (since angular momentum is T-odd), and 
is thus irrelevant in strong and electromagnetic matrix elements. The antisymmetric part is 
essentially a commutator of spin generators, and thus reduces to a single Pauli matrix; this 
is an example of an operator reduction rule. Thus, the only spin structure required up to 
the 3-body level is the operator with two Pauli matrices listed above. 

The quark charge operator can now either act upon one of the two quark lines on which 
the quadrupole spin tensor acts, or upon a third line. Using the general rule that the 
minimum number of gluon exchanges necessary to connect the quark lines of an n-body 
operator is n — 1, the former is a 2-body operator and thus has the suppression coefficient 
1/Nc, while the latter is a 3-body operator and thus is prefaced with a One finds that 
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only two operators appear up to the 3-body level: 



Q Nc C 



(3.1) 



where B and C are unknown coefficients of order unity, times a characteristic hadronic 
quadrupole size (in fm^). This implies in particular 4 all-orders relations among the 6 
quadrupole moments, as we see below. The coefficient powers of 1/Nc are the only explicit 
differences between this expression and that in the A^^^ = 3 case p9[] . 

In deriving expressions for the matrix elements of these two operators, it is of great 
advantage to note that the ground-state baryon representation — the analogue of the 56 
of SU(6) — is completely symmetric in combined spin and flavor indices. This means in 
particular that all quarks of a given flavor uoi d are completely symmetrized, and hence carry 
the maximal spin: = Nu/2 and Sd = Nd/2, where Ni is the number of quarks of flavor 
i in the baryon. Furthermore, one has the constraints + Nd = Nc and Nu — Nd = 
and for nonstrange states the rule I = J applies, where J = S„ + S(i is the total baryon 
spin. A complete list of compatible operators is thus = P, J3, I3, Sj^, and Sd'^- 

In the calculation one encounters operators such as Yli'" QiO'iz/2 = QuSuz + QdSdz- The 
values of Qu4 are fixed by the anomaly cancellation conditions of the standard model with 



gauge group SU{N^) x SU{2) x f/(l) to be |2 



Qu. 



c,t 



{N, + l)/2N„ Qd,s,b = {-N, + l)/2N,. 



(3.2) 



The calculation of all the necessary matrix elements for the 2-flavor case therefore requires 
only knowledge of two specific matrix elements for coupled angular momenta (ji, 72) J] 
of course, here ji and j2 stand for Su and Sd, the total spin angular momenta carried by the 
u and d quarks, respectively. These matrix elements are: 



(JM(ji js) \Jiz\ JM'{jij2)) = ^MSm'm 



1 + 



Jl(jl + 1)-J2(j2 + 1) 

J(J + 1) 



(3.3) 



and 



(JM(jU2)|Ji.| J-lM'(jU2)) 



M'M 



[J' - (jl - J2)][(jl + J2 + 1)2 - J']iJ' - M2) 



(2J + 1)(2J-1) 



(3.4) 



These expressions can be derived through elementary means. The first is a special case of 



the Wigner-Eckart theorem often called the projection theorem [33|. The second may be 



obtained solely through the manipulation of raising and lowering operators and is proved in, 
e.g., Ref. [Q. Since these theorems refer to the coupling of just two states of known angular 
momentum, they provide a complete answer only in the 2-flavor case. Of course, the 3-flavor 
case requires the coupling of an additional angular momentum (S^), which introduces 6j and 
9j symbols, as well as the careful consideration of additional symmetry properties of the 
baryon states. We do not pursue these topics here, but relegate their resolution to another 
paper |l20 |. 

The operators S^^^ are rank-1 tensors and therefore can connect states differing by up to 
one unit of total angular momentum J. Thus one obtains both diagonal and off-diagonal 
quadrupole transition matrix elements. Expressing Eq. (|3?I|) as Q = BOb + COq, one finds: 



{JJ\Ob\ J J) 
{JJ-1 \Ob\ J-lJ-1) 



+ 



1 



J 



{N, + 2){2Q 
J+1 



1) 4J 



3( J - 1) [(2g - 1)2 - 4 J2] [4J2 - (AT, + 2y 



2JNr 



{JJ\Oc\JJ) = -—(.J-- 



2J+1 
(iV, + 2)(2g-l) 

J + 1 



{JJ-l\Oc\J-lJ-l) 



3(J - 1) [{2Q - 1)2 - 4J2][4J2 - (iV, + 2)2] 



2 J 



(3.5) 



where Q is the total baryon electric charge. 



IV. RESULTS AND DISCUSSION 



Equations (|3.5|) evaluated for the quantum numbers of the 6 nonstrange states give the 
results listed in Table I. First note that these results agree with those in the = 3 case. 
Tables I and II of Ref. |]19[, except for our aforementioned coefficient factors of l/N^ and 



for B and C, respectively, and a factor of 2 accidentally neglected in the C terms 
of [|T9|. The matrix elements in were obtained by using explicit Nc = 3 spin-flavor 
baryon wave functions. 

Several features are immediately obvious. The diagonal matrix elements for Nc = 3 are 
just 4Q(35 + C)/9, while those in the ^ oo limit are 4/35/5 = 4(Q - 1/2)5/5, and 
the two off-diagonal elements are equal. These simple linear patterns arise from our model 
assumption that only one quark charge operator appears, and so the quadrupole operator 
transforms only as / = and 1. In particular, the 1 = 3/2 A's can be connected by operators 



8 



transforming as / = 0, 1, 2, or 3, and the 1 = 2 and 3 combinations must vanish for any A^^^ 
These combinations are: 



Qa++ - 3Qa+ + 3gAo -Qa- =0 (/ = 3), 
Qa++ - Qa+ - Qao + Qa- =0 (/ = 2). 



(4.1) 



Quadrupole transition operators between the / = 3/2 A and I = 1/2 N can transform as 
/ = 1 or 2, and since the latter does not appear in our model, it leads to the relation: 



Qa+p — Qaou, 



(4.2) 



for all Nc- The last of the 4 relations (6 degrees of freedom, 2 operators) is Nc dependent: 

Qao 2N^-3 \2{N^ + 5) 



5 Nr. 



Nr. 



(4.3) 



Qa+p 

and in particular vanishes for N^ = 3. 

In addition to combinations that vanish identically, one may consider combinations 
for which the leading in 1/Nc terms cancel. In particular, we note that generically, the 
quadrupole moments for large Nc are 0{N^), the sole leading contribution arising from the 
operator Ob- As discussed above, the diagonal matrix elements will prove extremely difficult 
to measure directly, and so we now seek to express all other quadrupole moments in terms 
of the one most easily measurable, Qa+p- In terms of the combination Q = 2\/2Qa+p/5, we 
find: 

1 



Qa++ 
Qa+ 
Qao 
Qa- 



+3Q 
+Q 
-Q 
-3Q 



1 + 
1 + 



N^ 
1 /39 

21 



19 C 
— + 20— 
2 B 



30 



.C 
B 



+ 



2N^ 

--(- 



10- 



B 



O 



lei 



(4.4) 



Phenomenological experience with the l/N^ expansion tells us that quantities for which 
corrections are only 0{1/N^) tend to agree well with their predicted central values. We 
should note, however, that the numerical coefficients of the 0{1/N^) terms in this case can 
be large for Nc = 3, depending upon the precise value of the ratio C/B. An interesting 



difference of prediction between this work and is the ratio Qa+IQa+p^ which here is 



2\/2lb + 0{1/Nl), and in Eq. (5) of [|T^] is -\/2; the latter prediction is obtained by setting 
C = (and of course Nc = 3). We defer detailed numerical analysis to such time as the 
3-flavor predictions are also in hand ^0 . 



Lastly, suppose that the baryon charge radii, which come from spin-spin terms in the 
baryon Hamiltonian, arise from the same source as the quadrupole operators; then the 
matrix elements for the two observables appear in fixed ratios. This occurs, for example, 
in one-gluon exchange picture, in which a multipole expansion of the baryon charge density 
operator p reads: 

Nc Nc 

P = AY^Qi- — '^Qi[2(Ti- (Tj - (3(Tiz(Tjz- (Ti- (Tj)] 

C 

~Jrl Yl Qk pCTj ■ (Tj - (ScTi^ - (Ti ■ (Tj)\ . (4.5) 



Note that the normalization of the spin-spin operator is just —2 times that used in Ref. [118 
Thus, the results obtained there can be carried directly over for comparison. In particular, 
the neutron charge radius is: 



c 



from which we immediately see that 



The factor on the r.h.s. is especially interesting since it equals 1 in both the A^^^ = 3 and 
Nc ^ oo cases (and in between never differs from unity by more than 1.2%). Thus, one 
expects that this relation should hold especially well. The A^^c = 3 version was first derived in 
Ref. using a constituent quark model. Using the value = — 0.113(3)(4) fm^ one 



predicts Qa+p = —0.0799(4) fm^, which agrees well with the value stated in Sec. II One 
could continue from here to predict all the baryon quadrupole moments using Eq. (|4.4| ). The 
reader should be reminded, however, that the assumption of one-gluon exchange is much 
more particular than the minimal assumption used to obtain Table I. 

V. CONCLUSIONS 

The 1/Nc expansion provides an additional handle on nonperturbative QCD phenomenol- 
ogy, through the observation that universes with odd Nc > 3 are very similar to our own. 
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Using a rigorous definition of constituent quarks obtainable in tlie large-A'^c limit and a 
minimal model ansatz — that quadrupole operators are proportional to the quark charge 
(single-photon exchange approximation) — we have shown that the 6 quadrupole moments of 
the nonstrange baryons can be described by just 2 distinct operators, leading to a number of 
relations that hold for all A^^^ and others that hold up to 0{1/N^) corrections. We have also 
seen that in a one-gluon exchange picture, all of these can be predicted using the measured 
value of the neutron charge radius r^. Detailed analysis of charge radii and quadrupole 
moments for the 3-flavor case is forthcoming. 
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TABLE I: Quadrupole moments of the nonstrange ground-state baryons for arbitrary Nc, = 3, 
and Nc oo. 



Qa++ 


B 


(+1) (iV2 + 2iV, + 5) + ^ (-f ) (iV2 - 8iVe + 5) 


§5 + 






Qa+ 


B 


(+§) (7V2 + 2iV, + 15) + j| (-1) (iV2 - 13iV, + 5) 




ic 




Qao 


B 


(-1) (iVe + 5)(Ar, - 3) + ^ (+1) (iVe + 5)(Ar, - 3) 









Qa- 


B 


(-1) (iV2 + 2iV, - 5) + ^ (+f ) (Ar2 - 3iVe - 5) 




-ic 




Qa+p 

QaOti 


( B o C A J {N,+5)(N,-l) 
(B o C \ /(iVc+5)(iVe-l) 


^{3B 


-2C) 
-2(7) 


71^ 
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